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The  performance  bounds  of  an  atmospheric  acoustic  array  operating  in  a 
turbulent  medium  with  fluctuations  described  by  a  von  K^rmdn  spectrum 
are  investigated.  This  treatment  considers  a  single  monochromatic  source 
and  a  line-of-sight  propagation  path.  The  primary  interests  are  in 
calculating  the  Cramer-Rao  lower  bounds  (CRLBs)  of  the  azimuthal  and 
zenith  angles  of  arrival  (AOAs)  and  in  observing  how  these  bounds 
change  with  the  introduction  of  additional  unknowns,  such  as  the 
normalized  propagation  distance  (to  wavelength),  turbulence  parameters, 
and  signal-to-noise  ratio  (SNR).  In  both  two  and  three  dimensions,  the 
CRLBs  of  the  AOAs  increase  significantly  for  large  values  of  the 
index-of-refraction  variance  and  normalized  propagation  distance.  For 
small  values  of  the  index-of-refraction  variance  and  normalized 
propagation  distance,  the  SNR  is  the  limiting  factor.  For  the 
two-dimensional  treatment,  the  estimate  of  the  AOA  wUl  decouple  from 
the  estimates  of  the  other  parameters  with  the  appropriate  choice  of  array 
geometry.  In  three  dimensions,  again  with  an  appropriate  choice  of  array 
geometry,  the  estimates  of  the  azimuth  and  zenitti  will  decouple  from  the 
estimates  of  the  other  parameters,  but  because  of  the  constraints  of  the 
turbulence  model,  they  will  remain  coupled  to  one  another. 
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1.  Introduction 


Acoustical  direction-finding  and  tracking  systems  will  likely  play  a 
prominent  role  on  the  future  battlefield,  where  situational  awareness  will 
be  a  key  factor  affecting  the  survivability  of  light-  and  medium-weight 
forces.  The  main  advantages  of  acoustical  sensors  are  low  cost,  small  size, 
passive  operation,  and  operational  capabilities  in  non-line-of-sight 
(non-LOS)  conditions.  However,  the  performance  of  acoustical  sensors 
does  strongly  depend  on  environmental  conditions. 

This  report  discusses  how  one  important  environmental  effect,  the 
scattering  of  soimd,  affects  the  ability  of  arrays  of  acoustical  sensors  to 
determine  the  bearing  angles  of  targets.  Scattering  occurs  when  sovmd 
interacts  with  turbulence  and  other  random  atmospheric  motions, 
creating  random  distortions  in  the  propagating  wavefronts.  As  the 
wavefronts  propagate  from  the  source  (target)  to  the  receiving  array,  they 
can  accumulate  substantial  random  variations  in  their  orientation  and 
intensity,  which  are  perceived  as  fluctuations  in  the  apparent  bearing 
angles  and  strength  of  the  source.  These  acoustic  phenomena  are 
analogous  to  scintillation  and  quivering  of  optical  images,  as  are  often 
observed  above  a  roadway  on  a  sunny  afternoon. 

The  study  of  the  scattering  of  waves  propagating  in  a  random  medium  is 
already  well  established.  Only  recenfiy,  however,  have  researchers  begun 
to  directly  incorporate  the  effects  of  the  scattering  into  performance 
predictions  for  acoustic  direction-finding  arrays.  Accounting  for  the 
strong  distortion  of  the  wavefront  by  atmospheric  turbulence  is  crucial  to 
obtaining  accurate  performance  characterization. 

The  performance  of  a  sensor  array  may  be  quantified  by  calculating  the 
mean  square  error  (MSE)  between  the  estimated  parameter  (such  as  the 
angle  of  arrival  (AO A)  or  propagation  distance)  and  its  actual  value.  The 
lower  bound  of  the  MSE  is  the  Cramer-Rao  lower  boimd  (CRLB),  which  is 
calculated  from  the  Fisher  information  (FI).  An  approach  that  directly 
incorporated  the  effects  of  a  random  medium  into  the  calculation  of  the  FI 
was  first  formulated  by  Song  and  Ritcey  [1]  to  calculate  the  performance 
boimds  of  AOA*  estimates  for  ocean  acoustics.  Using  the  general 

*In  this  report  we  make  a  distinction  between  the  AO  As  of  the  sound  wave  at  the  array 
and  the  actual  angles  of  bearing  (AOB)  of  the  source.  The  theoretical  treatment  in  this  re- 
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framework  of  Song  and  Ritcey  [1],  Wilson  [2]  calculated  the  performance 
bounds  on  AOA  estimates  for  atmospheric  acoustic  arrays. 

The  analysis  by  Song  and  Ritcey  [1]  considered  the  effects  of  a  random 
medium,  characterized  by  a  fluctuation  strength,  on  the  performance  of  a 
linear  acoustic  array.  The  total  received  signal  was  taken  to  be  that  from 
multiple,  uncorrelated,  plane-wave  source  signals  and  from  noise.  They 
assumed  that  both  the  signal  component  and  the  noise  component  had  a 
Gaussian  distribution  with  zero  mean.  The  covariance  matrix  for  the 
source  was  calculated  from  the  second  moment  of  the  sound  field  as 
predicted  by  the  theory  of  waves  propagating  in  a  random  medium.  For 
no  noise,  they  found  a  non-negligible  change  in  the  CRLB  of  the  AOA  for 
increasing  fluctuation  strength.  For  a  fixed  fluctuation  strength  and 
unknown  noise  parameter,  they  found  that  the  CRLB  of  the  AOA 
decreased  as  the  signal-to-noise  ratio  (SNR)  increased.  A  maximum 
likelihood  estimator  (MLE)  was  used  and  compared  to  the  CRLB:  the  MLE 
was  usually  well  within  20  percent  of  the  CRLB  for  five  independent  data 
sets.  The  authors  also  noted  that  the  estimator  calculated  from  a 
beamforming  method  with  the  use  of  a  Bartlett  processor  was  identical  to 
the  MLE  calculated  for  a  homogeneous  medium.  (Propagation  through  a 
homogeneous  medium  is  the  special  case  of  propagation  through  a 
random  medium  with  zero  fluctuation  strength.)  Comparisons  between 
the  two  estimation  methods  found  that  the  performance  of  the 
beamforming  method  was  considerably  worse  than  that  of  the  maximum 
likelihood  method  with  use  of  the  correct  fluctuation  strength.  Their 
findings  emphasize  the  necessity  of  correctly  incorporating  the  effects  of  a 
random  medium  on  the  wavefront. 

Wilson  [2]  used  the  same  general  approach  to  study  the  effects  of 
atmospheric  turbulence  on  the  performance  boimds  of  linear  acoustic 
arrays.  This  analysis  considered  a  single  monochromatic  plane-wave 
source  propagating  in  atmospheric  turbulence  and  a  LOS  propagation 
path.  The  CRLB  of  the  azimuth  was  calculated  with  the  assumption  that 
the  source  and  noise  components  of  the  received  signal  had  a  Gaussian 
distribution  with  zero  mean.  Gaussian,  von  Kdrm^,  and  Kolomogorov 
turbulence  models  were  examined.  It  was  found  that  the  array's 

port  applies  to  the  AOAs,  which  represent  the  orientation  of  the  wavefront  normal  when 
the  sound  reaches  the  array.  For  propagation  in  the  atmosphere,  the  average  horizontal  (az¬ 
imuthal)  AOA  is  usually  very  close  to  the  AOB,  thereby  making  acoustic  arrays  well  suited 
to  determining  the  horizontal  position  of  a  source.  The  situation  is  usually  quite  different 
for  the  vertical  (zenith)  AOA,  however.  In  this  case,  atmospheric  refraction  interferes  with 
the  ability  to  determine  the  AOB  from  the  AOA  by  bending  sound  waves  upward  or  down¬ 
ward.  Because  of  refraction,  most  existing  acoustic  systems  do  not  attempt  to  determine  the 
zenith  of  a  near-ground  source. 


performance  degraded  with  increasing  propagation  distance,  increasing 
frequency,  and  increasing  turbulence  strength. 

By  assuming  that  the  signal  of  interest  was  zero-mean  Gaussian,  both 
Song  and  Ritcey  [1]  and  Wilson  [2]  implicitly  treated  the  case  of  waves 
strongly  scattered  by  the  turbulence.  In  strong  scattering,  the  turbulence  is 
sufficiently  strong  and/or  the  wavefronts  propagate  sufficiently  far  that 
the  phase  of  the  received  signal  is  completely  randomized.  However,  for 
most  Army  tactical  scenarios  involving  acoustic  ground  sensors,  the 
propagation  distances  are  short  enough  that  the  signal  is  only  weakly 
scattered,  meaning  that  it  has  a  deterministic  mean  component.  This  report 
generalizes  the  treatment  of  the  previous  authors  by  providing  results 
valid  for  both  strong  and  weak  scattering. 

The  second  main  contribution  of  this  report  pertains  to  the  problem  of 
multiple  unknowns  in  the  FI.  In  Wilson  [2],  the  FI  was  calculated 
assuming  that  the  only  unknown  was  the  wavefront  AOA.  The 
source-receiver  propagation  distance  and  the  turbulence  parameters  were 
implicitly  assumed  to  be  known.  In  a  real  scenario,  this  information  may 
not  be  available.  In  this  report  we  extend  the  analysis  of  Wilson  [2].  The 
AOAs  (azimuth  and  zenith),  propagation  distance,  SNR,  and  turbulence 
parameters  are  treated  as  unlmowns.  The  first  and  second  moments  of  the 
scattered  sound  field  are  considered  for  a  turbulent  medium  with 
fluctuations  described  by  a  von  Karman  spectrum.  Both  a 
two-dimensional  (2D)  and  a  three-dimensional  (3D)  analysis  are 
performed  for  a  single  monochromatic  plane-wave  source.  A  LOS 
propagation  path  is  assumed.  The  primary  interests  of  this  report  are  to 
calculate  the  CRLBs  of  the  AOAs  and  to  observe  how  these  bounds 
change  with  the  introduction  of  additional  unknown  parameters. 
Therefore,  we  present  only  the  CRLBs  of  the  AOAs  in  this  report  and  leave 
the  discussion  of  the  CRLBs  of  the  other  parameters  for  a  later  report 
because  of  the  complexity  of  the  results. 

The  signal  processing  theory  is  described  in  section  2  and  the  von  Kerman 
turbulence  model  is  described  in  section  3.  The  fuU  theoretical  model  is 
derived  in  section  4.  In  sections  5  and  6,  the  CRLBs  of  the  estimates  are 
discussed  for  the  cases  of  no  turbulence  and  turbulence,  respectively. 
Results  are  given  in  section?  and  conclusions  are  drawn  m  sections.  A 
table  of  symbols  and  abbreviations  is  given  in  appendix  A. 
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2.  Signal  Processing  Theory 


2.1  Cramer-Rao  Lower  Bound 


We  wish  to  estimate  a  vector  parameter  ©  =  [©i  ©2  ...  © The 
minimum  MSE  of  an  unbiased  estimator  ©  about  its  actual  value  ©  for  a 
given  algorithm  may  be  calculated  from  the  Cramer-Rao  theorem  [3,4], 
which  gives 

((e„-e„)^)>  [j-'(e)]„^.  (D 

where  J (©)  is  the  Af  x  J\f  FI  matrix.  (An  estimator  is  said  to  be  unbiased  if 
and  only  if  ^0)  =  ©.)  The  FI  is  related  to  the  probability  likelihood 
p(x;  ©)  (the  probability  density  function  (PDF)  of  x  with  ©  as  a 
parameter)  by 

where  the  expectation  value  is  taken  with  respect  to  p(x;  ©)  and  the 
derivatives  are  evaluated  at  the  true  value  of  ©.  For  some  function  the 
estimator  ©  =  ^(x)  can  attain  the  CRLB  if  and  only  if  [3] 


91np(x;  ©) 

a© 


=  J(©)[^(x)-©]  . 


(3) 


The  likelihood  function  for  real  parameters  of  a  complex  Gaussian  PDF 
with  covariance  matrix  Cx(©)  and  mean  /x(©)  may  be  written  [3] 

°  d^t [c,(e)]  {"  ~ >‘<^>1'  c;'(e)ix-,,(e)]}.  (4) 

Its  corresponding  FI  is 


=  tr 


d/jL^ 


c-i 


(5) 


where  the  functional  dependence  has  been  dropped.  If  there  are  M 
independent  and  identically  distributed  data  sets,  the  likelihood  function 
is  the  product  of  M  identical  distribution  functions,  and  hence  the  FI  is  M 
times  the  quanitity  given  in  equation  (5).  Let  there  be  n  elements  in  the 
sensor  array.  Then  C  is  an  n  x  n  matrix,  and  is  a  column  vector  of  length 
n.  Let  the  subscripts  A,  i/  €  [1,  ... ,  Af]  be  the  indices  on  the  parameters 
and  i,  j  e  [1,  . . . ,  n]  be  the  indices  for  the  elements  of  the  sensor  array. 
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Let  us  define  ^^7  sometimes  refer  to  the  square  root 

of  the  CRLB  as  SQRTCRLB.)  The  minimum  value  of  al  is  al^  =  l/  J^,  i-e., 
the  CRLB  when  ©j,  is  tiie  only  unknown.  As  the  number  of  unknowns 
increases,  al  will  increase.  This  is  readily  seen  by  examining  the  FI  for  two 
unknowns.  Suppose 

J  = 


Jll  Ji2 
J21  J22  J 


(6) 


Then,  as  the  FI  is  symmetric  for  real  parameters. 


a!  =  (Jll  -  Jf2/J22) 


-1 


(7) 


Only  if  J12  =  0  do  (tJ  2  =  <^10  2o'  estimates  of  ©i  and  ©2  are  said  to 

be  uncoupled.  The  quantity 

^2  J2 

C  =  l-  ^  =  l-  -%  =  (8) 

(jf  (J2  J11J22 

provides  a  measure  of  the  strength  of  the  coupling  between  the  estimates 
of  ©1  and  ©2:  if  C  =  0,  the  estimates  are  imcoupled;  if  C  <  1/  the  estimates 
of  ©1  and  ©2  are  nominally  coupled;  and  if  =  1,  there  is  not  enough 
information  to  obtain  an  estimate  of  either  ©1  or  ©2  and  the  CRLB  of  each 
is  infinite. 


2.2  Signal  Model 

Consider  an  acoustic  array  with  n  sensors.  We  assume  that  the  signal  at 
each  sensor  results  from:  (1)  the  wave  that  has  propagated  from  the  source 
of  interest,  with  ip  and  6  as  the  azimuthal  and  zenith  AOAs,  and  (2) 
random  noise.  Let  p{p,6,  t)  and  N (t)  be  the  time-varying  complex 
envelopes  of  the  two  contributions,  respectively.  These  column  vectors 
have  n  elements,  one  element  corresponding  to  each  sensor.  The  source 
contribution  is  time  dependent  because  of  the  random  tiubulent  effects. 
The  noise,  which  is  also  time  dependent,  may  result  from  wind  noise  or 
other  competing  acoustic  sources.  The  total  received  signal  is  thus 

s{ip,B,t)  =  p{(p,e,t)  +  N{t).  (9) 

Let  us  assume  that  the  source  signal  and  the  noise  are  uncorrelated.  Then 
if  the  noise  signals  at  the  sensors  are  mutually  vmcorrelated  and  have  an 
equal  variance  of  cr^  ,  the  total  covariance  is  the  sum 

C  =  Cp  +  C^  =  Cp  +  crlln,  (10) 
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where  the  subscript  p  refers  to  the  pressure  field  of  the  sound  wave  of 
interest  and  N  the  noise.  We  assume  that  the  noise  has  a  Gaussian 
distribution  with  zero  mean.  Now  by  definition 

=  {{Pi-  Pi){Pj  -  Pj)*  )  ^  (piP*)  -  pip*,  Pi  ={  Pi),  (11) 

where  (•)  indicates  the  ensemble  average.  We  take  the  values  of  {piP*  ) 
and  Pi  from  results  in  the  literature  for  acoustic  wave  propagation  in  a 
moving  random  medium  as  discussed  in  section  3.1. 


3.  Turbulence  Model 


Ideally,  we  should  calculate  the  FI  from  the  product  of  the  probability 
Hkelihoods  of  the  pressure  field  p  and  the  noise  N.  However,  exact 
solutions  for  p  do  not  exist  (nor  for  p(p)),  but  solutions  to  the  first  and 
second  moments  do  exist. 

The  rigorous  solution  for  the  pressure  field  associated  with  a  sound  wave 
propagating  through  a  moving  medium  can  be  obtained  as  a  closed  set  of 
fluid  dynamic  equations  [5].  A  pair  of  approximations  are  commonly 
applied  to  this  equation  set  to  make  the  problem  tractable.  The  parabolic 
approximation,  valid  for  small-angle  propagation,  provides  reduction  to  a 
single  equation.  The  Markov  approximation,  which  assumes  that  the 
turbulence  field  has  vanishing  correlation  in  the  propagation  direction, 
subsequently  allows  the  statistical  moments  of  the  sound-pressure  field  to 
be  obtained  in  a  closed  form.  These  parabolic  and  Markov  approximations 
are  valid  in  the  far  field  of  the  source,  at  observation  points  near  the 
propagation  axis  <  1  where  ^  is  the  angle  off  the  propagation  axis), 
and  for  £  A  >  ^,  where  A  is  the  wavelength  and  C  and  (.  are  the  outer 
(integral)  and  iimer  length  scales  of  the  turbulence,  respectively.  We  use 
the  solutions  for  the  first  and  second  moments  of  the  sound  pressure  as 
given  by  Ostashev  [5].  These  results  are  outlined  in  the  following  section. 

3.1  Solution  for  the  First  and  Second  Moments 

Suppose  that  a  sound  wave  is  propagating  along  the  x-axis.  The  parabolic 
and  Markov  approximations  predict  that  the  moments  of  the  soimd 
pressure  will  attenuate  in  a  simple  exponential  manner.  Defining  the 
observation  point  as  r  =  [x,  y,  z]  and  the  acoustic  wave  number  as  k, 

|k|  =  A:  =  27r/A,  the  equation  for  the  first  moment  is  then 

<p(r)>  =  e-^"pH(r),  (12) 

where  7  is  the  extinction  coefficient  for  the  first  moment  and  Ph  (r)  is  the 
soimd  field  in  the  absence  of  random  inhomogeneities.  For  an  incident 
plane  wave,  pH(r)  =  po  =  po  where  po  is  a  real-valued  constant. 

We  also  define  a  saturation  parameter  O,  as 

D  =  1  -  .  (13) 
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The  first  moment  represents  the  unscattered  (deterministic)  part  of  the 
wave  field.  When  <  1,  the  saturation  parameter  is  close  to  0,  scattering 
is  small,  and  the  signal  at  a  single  sensor  exhibits  little  variability.  When 
7a;  »  1,  1,  and  the  scattered  part  of  the  field  dominates. 

We  wish  to  know  the  second  moment  for  two  points  near  the  x-axis, 
designated  ri  =  [x,  1/1,  zi]  and  r2  =  [x,  y2,  Z2].  For  a  fixed  sensor 
separation,  the  second  moment  also  undergoes  an  exponential 
diminishment: 


(  p  (ri)  p*  (r2)  >  =  pn  (ri)  p*  (r2)  e  ,  (14) 

where  p  =  ri  -  r2  is  the  sensor  separation  vector  (transverse  to  the 
propagation  direction)  and  a  is  the  extinction  coefficient  for  the  second 
moment. 

The  extinction  coefficients  depend  on  the  structure  of  the  random 
medium.  The  result  for  the  second  moment  is  [5-8] 

a(p)  =  27rk'^  [/(O)  -  /(p)]  (15) 

in  which  /  is  the  two-dimensional  (or  projected)  correlation  hmction  for 
the  sound-speed  fluctuations.  The  extinction  coefficients  for  the  first  and 
second  moments  are  related  by 

7  =  a(oo)/2 .  (16) 

For  most  random  media,  including  turbulence,  Q;(p)  initially  increases 
monotonically  with  increasing  p,  but  when  p  exceeds  C,  a{p) 
asymptotically  approaches  a  constant  value  [6].  Since  /(p)  0  in  the  limit 

p  — >  00,  this  constant  value  is  simply  27,  given  by 

27  =  27r/cV(0)=2c2fe2£,  (17) 

where  is  the  index-of-refraction  variance.  Hence  the  second  moment 
will  initially  decrease  with  increasing  p,  but  will  eventually  "saturate"  at  a 
fixed  minimum  value. 

3.2  von  Karmin  Turbulence  Model 

Typical  acoustic  sensor  arrays  used  by  the  Army  have  a  sensor  spacing 
larger  than  the  height  of  the  array  from  ground.  As  such,  performance  of 
these  arrays  is  affected  by  the  large  eddies  of  the  energy-containing  (or 
source)  subranges  of  the  turbulence  spectrum.  This  dependence  on  the 
large  eddies  distinguishes  typical  acoustical  systems  from  optical  ones. 
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whose  performance  depends  primarily  upon  the  smaller  scale  eddies  in 
the  inertial  and  dissipation  subranges.  In  this  report,  we  adopt  the  isotropic, 
homogeneous  von  K^rm^  turbulence  model,  since  it  is  relatively  simple 
and  behaves  fairly  realistically  in  the  energy-containing  subrange.  And, 
unlike  the  more  commonly  used  Gaussian  models,  the  von  Kdrmdn  model 
accurately  describes  the  inertial  subrange  of  the  turbulence  spectrum. 

The  von  Karman  form  for  the  two-dimensional  correlation  function 
actually  depends  on  whether  the  fluctuations  in  soimd  speed  are  induced 
by  a  scalar  field  (temperature  or  humidity)  or  a  vector  field  (wind 
velocity).  The  fimctions  for  a  scalar  field  /«  and  a  vector  field  fv  are  [6] 


fs{p, 

fv^Pi^  jO 


2<;H 

2<;H 


(18) 

(19) 


where  Z  is  a  characteristic  length  scale  (equal  to  r(l/3)  r  ^(5/6)  tt  C), 
r(x)  is  the  gamma  function,  and  K^^ix)  is  the  modified  Bessel  function  of 
order  v. 
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4.  Full  Theoretical  Model 


We  now  wish  to  apply  the  results  presented  in  the  previous  sections  to 
calculate  the  performance  boimds  of  an  acoustic  sensor  array  operating  in 
a  turbulent  medium  with  fluctuations  described  by  a  von  Karmdn 
spectrum.  To  do  so,  several  approximations  must  be  made.  First,  since  a 
closed  form  expression  for  p  (p)  does  not  exist  in  the  solution  to  the 
parabolic  equation,  we  make  the  approximation  that  p  (p)  has  a  Gaussian 
distribution  whose  covariance  matrix  is  described  by  the  first  and  second 
moments  of  p  as  presented  in  the  section  3.1.  Second,  the  solution  of  the 
second  moment  given  in  section  3  is,  strictly  speaking,  only  valid  in  the  far 
field  for  two  observation  points  that  have  an  equal  distance  parallel  to  the 
wavefront  normal  n  and  have  small  components  transverse  to  fi.  This  is 
an  unhkely  scenerario,  however,  when  one  performs  source  localization 
with  an  acoustic  array. 

Let  r  be  the  vector  from  the  center  of  the  array  normal  to  the  plane  of  the 
source  (parallel  to  n),  so  that  r  is  the  propagation  distance  of  the 
wavefront  to  the  center  of  the  array.  Let  =  [x[,  y\,  z']  be  the  vector  from 
the  center  of  the  array  to  the  ith  sensor.  We  define  Rj  =  r  +  r^,  so  that  i?|  is 
the  propagation  distance  of  the  wavefront  to  the  ith  sensor.  Let  be  the 
separation  between  the  ith  and  jth  sensors,  =  R,  -  -  r' ,  and 

let  pf-  be  the  separation  transverse  to  n.  See  figure  1.  Strictly  speaking 
equation  (14)  is  valid  only  for  p^  =  pjj.  However,  we  approximate 

We  have  two  approximations  here.  First,  we  assume  that  the  smaU  angle 
approximation  holds,  so  that  pjj  ^  p^j  and  hence  a{pjj)  ~  a(pij).  Second, 
we  take  x  in  equation  (14)  to  be  r  =  |r|.  This  is  valid  provided  that  for 
every  i  and  j, 

e-°{pij)r  ^  ^-a{Pij)Ri  ^  ^-a{p^J)R'j  ^2^) 

Caution  must  be  exercised  when  applying  these  approximations.  For 
example,  suppose  that  there  are  four  sensors  arranged  in  a  rectangular 
configuration.  If  n  is  in  the  plane  of  the  sensor  array,  then  there  exist  an  i 
and  j  such  that  pij  pj-y  Whereas,  if  n  is  nearly  perpendicular  to  the  array 
plane,  for  every  i  and  j,  pij  «  See  figure  2.  Therefore,  when 
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Figure  1.  Coordinate  system  depicting  vector  from  array  center  (origin)  to  ith  sensor  r^, 
between  ith  and  jth  sensors  =  r^  —  r' ,  from  array  center  normal  to  source  plane  r  (so 
that  r  is  propagation  distance  to  array  center),  and  R,  =  r  +  r^.  Closed  circles  represent 
sensors  and  open  circle  represents  point  at  which  r  is  normal  to  source  plane. 


considering  problems  strictly  in  two  dimensions,  we  use  a  linear  array. 

The  use  of  other  array  geometries  requires  a  different  theoretical 
formulation  than  presented  here.  Similarly,  when  considering  problems  in 
three  dimensions,  we  avoid  the  use  of  three-dimensional  cirray 
configurations  such  cubes. 

4.1  Covariance  Matrix 

Based  on  the  assumptions  just  presented,  we  take  the  first  and  second 
moments  to  be 

fii  =  (22) 

(23) 

where 

(24) 

(25) 

(26) 

and  where  r,  R^,  and  Pq  are  as  previously  defined*.  The  phase  of  the 
source  %  is  an  unknown  parameter  and  must  be  included  in  this  treatment 

‘Note  that  we  use  the  convention  that  fii  is  the  ith  component  of  the  vector  /x  (i.e.,  the 
value  of  the  first  moment  of  the  sound  field  at  the  ith  sensor)  and  is  given  by  p,  =  .  In 

this  manner,  |pi|  =  p  and  |p|  =  p  =  ^/njl. 


ft  =  Poe 

=  k  •  R,  +  X 
4>ij  =  (fi  (fj  ) 
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(a)  Two  dimensions 


(b)  Three  dimensions 


Figure  3.  Coordinate  system  in  (a)  two  dimensions  and  (b)  three  dimensions. 


as  we  are  assuming  a  nonzero  mean.  We  thus  write  the  elements  of  the 
total  covariance  matrix  (equation  (10))  as 


2  ^  2  I  2 

Cii  =  Po  —  fJ-  + 


Cij  — 


p2g  a(pij)r_p2 


(27) 

(28) 


Let  us  define 


r  =  r  [cos  sin  0,  siny>sin0,  cos0] 

i-i  =  [x'i,  y'i,  4]  > 


(29) 

(30) 


so  that  ip  is  the  azimuth,  0  is  the  zenith  (or  declination),  and  9'  =  tt/2  -  Ois 
the  elevation.  See  figure  3.  Thus 

k  •  Ri  =  fcf  •  (r  +  r[)  =  A;  (r  +  cos  9?  sin  0  +  sin  sin  0  +  z'  cos  9)  .  (31) 


The  FI  may  now  be  readily  calculated  for  those  parameters  we  wish  to 
consider  as  unknowns.  While  the  derivatives  of  the  covariance  matrix  and 
mean  are  straightforward,  the  derivatives  with  respect  to  the  turbulence 
parameter  I  are  tedious.  For  the  interested  reader,  these  are  presented  in 
appendices  B  and  C. 


4.2  Numerical  Issues 

Numerical  difficulties  may  arise  in  trying  to  invert  J.  Consider  the 
diagonal  terms  of  ihe  FI  matrix.  The  contributions  from  the  propagation 
distance  r  and  characteristic  length  scale  I  may  differ  by  more  than  10 
orders  of  magnitude.  This  scale  discrepancy  may  also  occur  for  other 
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parameters.  We  construct  a  diagonal  matrix  D  whose  elements  are  the 
square  root  of  the  diagonal  elements  of  the  FI  matrix,  Dxu  =  The 

matrix  D  JD  is  now  invertible.  The  CRLB  of  the  unknown  parameters  may 
then  be  determined  from 


(  (0^  -  0^)2  )  =  [d  (DJD)"^  D 


(32) 


The  efficiency  of  the  estimator  is  maintained  by  this  process,  as  it  is 
equivalent  to  a  linear  transformation  of  0.  (Linear  transformations 
preserve  the  efficiency  of  an  estimator  [3].) 
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5.  Discussion:  No  Turbulence 


We  can  gain  insight  into  the  estimation  of  the  AOAs  by  first  examining  the 
limiting  case  of  no  turbulence.  In  the  absence  of  turbulence,  we  simply 


have 


C  =  cr^I„  and  /*  =  po 


(33) 


For  parameters  ©a  and  Qj,  of  which  cr^  is  independent,  the  Fisher 
information  is  determined  from  the  second  term  in  equation  (5)  and  is 


given  by 


2M  dfj, 
^  d0A  50, 


-  oafPo  ^ 


5.1  Angle-of-Arrival  Estimate  in  Two  Dimensions 

In  two  dimensions,  4>i  =  kr  +  kx\  cos  tp  +  ky[  sin  <p  +  x-  assume  r  is 
known,  if  not,  we  can  make  the  substitution  x  — >  +  X  below,  and 

the  results  will  stiU  hold.  We  have 


(35) 

(36) 

(37) 


The  estimates  of  cp  and  x  will  decouple  when  =  0  X)  2/i  = 

i.e.,  if  the  origin  is  taken  to  be  at  the  center  of  the  array. 

Note  that  none  of  the  elements  of  the  FI  is  dependent  upon  the  value  of  x- 
This  is  expected  as  it  is  the  arbitrary  phase  of  the  source,  and  its 
dependence  should  eventually  cancel.  Close  inspection  of  the  second  term 
of  equation  (5)  reveals  that  its  dependence  in  the  FI  should  cancel 
regardless  of  whether  we  consider  turbulence  or  not. 
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5.2  Inclusion  of  x  for  Non-Zero  Mean  Case 

It  was  noted  earlier  that,  as  we  are  considering  the  case  with  a  non-zero 
mean,  the  phase  of  the  source  x  in  equation  (25)  must  be  included.  The 
reason  is  most  readily  seen  by  examining  the  no  turbulence  case.  For 
simplicity  let  us  consider  a  line  array  along  the  y-axis  and  let  us  make  no 
assumptions  about  the  choice  of  the  coordinate  system.  First  suppose  that 
<f)i-kr  +  ky[  sin  ip  where  ip  is  unknown  and  r  is  known.  Then 


=  2Mk^  cos^ 


(38) 


As  —  1/ we  see  that  we  could  choose  a  coordinate  frame  infinitely 
far  away  from  the  sensor  array  so  that  the  CRLB  of  y?  is  zero. 

As  this  is  illogical,  let  us  now  suppose  that  <f>i  =  kr  +  ky[  sin  y>  -h  x  where 
both  ip  and  x  are  unknown  and  r  is  known.  The  CRLB  of  p  is  now 


^  t2  \  ^  2 

Jipx\ 

,  '^XX 


—  „ji 

pI 


pI 


Mk^ 


cos 


r 

”  1 

/  n  \  2 

> 

J  2M k^  cos^  p 

i2  1 

-n 

l~l 

fe')j 

> 

V. 

1-^ 

j 

n 


hO 


(39) 


Now  is  only  dependent  upon  the  separation  of  the  sensors,  regardless 
of  the  choice  of  coordinate  system. 

Second,  let  us  consider  that  (f)^  =  kr  +  ky[  sin  p  where  r  and  p  are  both 
unknown.  Then  is  given  by  equation  (38)  above  and 


=  2Mnk‘^^ 

(T^ 

(40) 

N  \ 

N 

- 

^  dr  dp 

=  2Mk^co^p^f^y\. 

(41) 

Let  us  choose  the  origin  to  be  at  the  center  of  the  array  so  that  the  estimates 
of  r  and  p  are  imcoupled.  Then  =  {2Mnk'^pl/al)~^ ,  i.e.,  we  have  the 
unplausible  result  that  the  propagation  distance  can  be  determined  to 
within  the  noise-to-signal  ratio  (NSR)  for  a  plane  wave  with  a  single 
sensor  array.  Clearly,  this  bound  could  never  be  attained  in  practice. 

Let  us  thus  suppose  that  =  A:r  +  ky[  sin  y;  +  x  where  r,  p,  and  x  are  all 
unknown.  Again  let  us  chose  the  origin  to  be  at  the  array  center,  so  that 
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the  estimate  of  ip  is  uncoupled  from  the  estimates  of  r  and  x-  The  FI  matrix 
is  then  block  diagonal,  and  the  CRLBs  of  r  and  x  can  be  determined  from 
the  submatrix 

J  =  2Mn^ 

N 

where  here  Jn  corresponds  to  Jrx-  Because  this  matrix  is  singular,  we 
have  the  expected  result  that  the  propagation  distance  and  phase  cannot 
be  determined  simultaneously  for  a  plane  wave. 

Clearly  from  both  examples  presented,  the  phase  angle  x  must  be 
included  in  (pi  for  the  non-zero  mean  case. 

5.3  Angle-of-Arrival  Estimates  in  Three  Dimensions 

Now  let  us  consider  the  CRLB  of  the  azimuth  and  zenith  in  three 
dimensions.  (We  note  that  we  may  interchangeably  refer  to  the  CRLB  of 
the  zenith  {$)  and  elevation  (O'),  for  as  they  are  related  by  a  linear 
transformation,  their  CRLBs  are  the  same.)  Suppose  that 
^i  =  kr  +  kx'i  cos  (p  sin  6>  -h  ky'^  sin  ^  sin  6>  -h  kz'^  cos  6>  -h  X  where  ip,  0,  and  x 
are  unknown  and  r  is  known.  We  have 


k'^  k 
k  1 


(42) 


2  ^ 

=  2Mk^  sin^  0^  ^  (-a:  ■  sin  p  +  y'i  cos  pY 

i=l 

2  ^ 

Jqq  =  2Mk^^  ^2  {^i  pcos0  +  y'i  sin pcosO-  z'^  sin O)" 


N  i=l 


Jyy  =  2Mn^ 
^2 


(43) 

(44) 

(45) 


2  ''' 

J^e  =  2Mk^  sin  6^  ^  (yf  -  ^  cos  (f  sin  (fcos6  +  x[y[  (cos^  -  sin^  (p)  cos  0 


N  i=l 


+  x[z[  sin  9?  sin  0  —  y[z[  cos  p  sin  0 


2  ^ 

=  2Mksm6^  ^  {x'i  sin p  +  y'i  cos p) 


N  i=l 


9  n 

Po 


Jsy.  =  2Mk^  ^  {x'i  cos  p  cos  0  +  y'i  sin  p  cos  0  —  z'^  sin  6) 


(46) 

(47) 

(48) 


N  i=l 


(Again  note  that  the  FI  is  independent  of  the  value  of  x-)  For  the  estimate 
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(49) 


of  X  to  decouple  from  the  estimates  of  ip  and  0 

n  n  n 

i=l  i=l  t=l 

and  for  the  estimates  of  ip  and  9  to  decouple  from  one  another 

=  and  =  0  .  (50) 

i=l  i=l  i=l  i=i  i~i 

Symmetric  array  configurations  such  as  a  circular  array  with  sensors 
placed  at  equal  angular  intervals  or  a  rectangular  grid  with  sensors  placed 
at  the  lattice  points  meet  these  requirements  provided  that  the  origin  is 
taken  to  be  at  the  center  of  the  array  and  that  the  array  is  located  in  the 
xy-plane.  (We  should  note  that  Nielsen  [9]  has  performed  a  similar 
analysis  for  a  multiple-frequency,  far-held,  sine  wave  signal  imbedded  in 
Gaussian  noise.  The  conditions  he  foxind  for  the  estimates  of  all  the  angles 
to  decouple  are  the  same  as  for  the  simple  case  presented  here.  Among  the 
literature  which  examine  array  conhgurations  that  result  in  the 
decoupling  of  the  angle  estimates,  [10-12]  may  be  of  interest  to  the  reader.) 

If  conditions  (49)  and  (50)  are  met,  then 

=  2Mk‘^  sin^  0^  ^  xf  (51) 

'^N  i=l 

J00  =  2Mk^^  cos^  6  +  sin^  0^  .  (52) 

N  i=l 

Note  that  the  CRLB  of  ip  is  singular  at  0  =  0.  Clearly  this  choice  of 
coordinate  system  is  benehcial  for  estimations  of  angles  that  are  not  near 
the  z-axis.  Recall  when  we  introduded  the  turbulence  model  that,  as  we 
are  using  the  small  angle  approximation,  we  are  interested  in  angles  near 
the  propapation  axis.  We  therefore  will  not  be  able  to  uncouple  the 
estimates  of  the  zenith  and  azimuth  when  considering  turbiilence,  as  we 
will  have  to  choose  the  array  plane  to  be  either  the  xz-  or  yz-plane. 

In  the  following  sections  when  we  consider  the  case  with  turbulence,  we 
will  take  the  array  to  be  in  the  yz-plane  and  the  wavefront  to  propagate 
near  the  x-axis.  Let  us  therefore  examine  the  case  of  no  turbulence  for  this 
configuration.  We  assume  that  equation  (49)  holds,  so  that  the  estimates  of 
p  and  0  are  lancoupled  from  the  estimate  of  x-  Moreover,  consider  an  array 
geometry  such  that 

X^y'z'  =  0  and  =  =  Y'\  (53) 

i=l  i=l  i=l 
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(a  square  grid  is  one  example).  Then 

J^if,  =  2Mk^Y'^^  sin^  6  cos^  tp  (54) 

Jee  =  2Mk^Y'‘^^  (sin^  (p  cos^  6  +  sin^  9)  (55) 

2 

J^e  =  2Mk‘^Y'^  ^  cos  p  sin  p  cos  0  sin  6  .  (56) 


Note  that  the  coupling  between 
dependent, 

72 

: 

JipipJ$0 


the  estimates  of  p  and  9  is  only  angularly 


sin^  p  cos2  9 
sin^  9  +  sin^  p  cos^  9 


(57) 


It  follows  that 


2Mk^Y'^pl  sin^  9 

2  <7^  sin^  9  +  sin^  p  cos^  9 

2MA;2y'2p2  (.Qg2  ^p  gjjj4  Q 


(58) 

(59) 


While  al  is  only  dependent  upon  9,  is  dependent  upon  both  p  and  9. 

As  we  are  considering  propagation  near  the  a;-axis,  it  is  useful  to  write  the 
results  in  terms  of  the  elevation  (9')  as  opposed  to  the  zenith  (0).  In  this 
marmer,  we  may  use  the  small  angle  approximation  to  the  power  series 
expansion  in  both  p  and  9'.  The  coupling  coefficient  is 

^  sin2(^siii2  9' 
cos^  0'  +  sin^  ip  sin^  0' 

»  .  (61) 


Note  that  for  either  97  =  0  or  =  0,  C  =  0  and  the  estimates  of  p  and  9'  are 
imcoupled  as  expected.  We  also  see  that  for  small  angles  the  coupling 
between  the  estimates  of  p  and  9'  is  weak.  In  fact,  at  ^  =  0'  =  15°, 

=  4.79  X  10~^.  The  CRLBs  of  9'  and  p  are 


crz 


2Mk^Y'^pl 


sec 


1  +  + 


(62) 

(63) 
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and 


(64) 


< _ ( 

V 


l  +  0'Vy^2  + V+?¥’'  +  20'V 


2Mk‘^Y'^pl  V'  '  "  '  3'  ’3’ 


(65) 


Note  that  (t^{p,  O')  =  <^^{0',  (p)  holds  up  to  Ath  order  in  the  expansion;  for 
V’  =  =  15°,  there  is  a  0.13  percent  difference  between  the  exact 

expression  and  the  Ath-order  approximation,  and  a  1.5  percent  difference 
for  the  2nd-order  approximation. 


6.  Discussion:  Turbulence 


6.1  Angle-of-Arrival  Estimate  in  Two  Dimensions 

Let  us  consider  the  CRLB  of  the  AOA  ip  in  two  dimensions  when 
incorporating  turbulence.  In  the  presence  of  turbulence,  we  find 
(numerically)  that  choosing  the  origin  to  be  at  the  center  of  the  array  is  no 
longer  a  sufficient  condition  for  the  estimates  of  the  azimuth  and  phase 
angle  to  decouple.  In  addition,  we  must  have  a  symmetric  array 
configuration  such  as  a  linear  array  with  mirror  symmetry  about  the 
origin,  (note  that  the  stronger  condition  of  uniform  spacing  is  not  needed). 
The  estimates  of  cp  and  x  will  also  be  xmcoupled  for  a  set  of  parallel  linear 
arrays  with  imiform  sensor  spacing,  however,  the  condition  for  the  small 
angle  approximation,  given  in  equation  (20),  is  violated  with  this  array 
geometry.  The  results  are  independent  of  the  value  of  %/  as  we  noted 
earlier  they  should  be. 

In  Wilson  [2],  the  CRLB  of  the  AOA  was  analysed  in-depth  for  the 
two-dimensional  case  assuming  a  zero  mean  and  a  single  unknown 
parameter  ip.  We  now  wish  to  examine  the  difference  when  the  mean  is 
incorporated.  For  simplicity,  let  us  begin  by  examing  a  simple 
two-element  array  and  considering  the  unknowns  p  and  %.  We  assume 
that  all  other  parameters  are  known.  In  appendix  D,  the  FI  for  a  two-point 
sensor  array  is  derived.  In  deriving  equations  (D-l)-(D-7),  no  assumptions 
are  made  about  the  form  of  the  derivatives  of  the  covariance  matrix  or  the 
mean  vector.  In  deriving  equations  (D-8)  and  (D-9),  it  is  assumed  that 
da/dQx  =  djS/dQx  =  dp/dOx  =  0,  where  a  =  Cn  =  C22  =  Po  "  in 

equation  (27)  and  =  IC12I  =  |C'2i|  =  in  equation  (28).  As 

noted  before,  the  estimates  of  the  two  parameters  ©a  and  0„  will  be 
uncoupled  if  Jxu  =  0. 

Suppose  that  Qx  =  P  and  0„  =  x-  By  substituting  the  derivatives 
presented  in  appendix  B  into  equation  (D-9),  we  find 

^  _  2M‘^p?  (d(i)i  ,  d(f)2\  _  [-  {x[  -|-  X2)  sin p  +  {y[  -I-  j/2) cos p] 

a  +  /3  {dp^  dp)  ~  a  +  /3 

Therefore,  the  estimates  of  p  and  %  will  be  tmcoupled  provided  =  -Xj 
and  y[  -  -y'^.  Let  x[  =  xg  =  0  and  y[  =  -yi^  =  d/2.  Then  the  CRLB  of  p  is 
which  from  equation  (D-8) ,  is 
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CM  ^ 


(a)//>.  =  5  (b)//7.  =  100 


Index-ol-refraclion  variance  c‘  Index-of-ref'raction  variance  c" 


Figure  4.  Percent  difference  of  cr,^  calculated  with  a  nonzero  mean  and  with  a  zero  mean 
versus  is  in  (a)  and  (b)  for  different  values  of  1/ A.  Corresponding  saturation  parameter  is 
in  lower  figures  (c)  and  (d).  Calculations  are  in  2D  for  (^  =  0,  r/A  =  250,  d/A  =  0.5,  n  =  5, 
and  a  scalar  von  Kdrman  spectrum. 


_ _ 

Mk'^d?  cos^  ip  [2/3^  -p  ^2 

pt  [1  -  -  2plt^  [1  -  +  al  (2pg  -  +  gg.) 

Mk'^cP  cos^  +  p'^p^  [l  —  +  p?a‘^  } 

For  p  =  0we  recover  the  result  presented  in  Wilson  [2].  Whether  or  not 
the  mean  is  incorporated,  when  ip  =  0,  the  source  is  broadside  to  the 
sensors,  and  the  performance  is  maximal.  When  p  =  7:12,  the  source  is  in 
line  with  the  sensors,  and  the  solution  is  singular.  Even  though  y;  =  7r/2  is 
not  consistent  with  the  small  angle  approximation,  the  solution  should  be 
singular  aitp  =  7r/2. 

In  figxires  4  to  6  the  percent  difference  between  with  and  without  the 


(67) 

(68) 
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(b)  ?2=10-' 


(a)  =  1 0  ^ 


(c)  <;-=10-" 


(d)  q-=io-' 


Figure  5.  Percent  difference  of  a.^  calculated  with  a  nonzero  mean  and  with  a  zero  mean 
versus  Z/ A  is  in  (a)  and  (h)  for  different  values  of  Corresponding  saturation  parameter  is 
in  lower  figures  (c)  and  (d).  Calculations  are  in  2D  for  <p  =  0,  r/A  =  250 ,  d/A  =  0.5 ,  n  =  5, 
and  a  scalar  von  Kdrmdn  spectrum. 


mean  term  is  presented  using  a  scalar  von  Karm^  spectrum.  A  linear 
array  with  5  sensors  spaced  in  equal  intervals  of  d  is  considered.  The 
upper  graphs  (a)  and  (b)  in  figure  4  depict  fi\e  percent  difference  versus 
the  index-of-refraction  variance  for  fixed  values  of  r/ A,  l/X,  and  d/ A 

and  two  values  of  the  SNR.  The  lower  graphs  (c)  and  (d)  give  the 
corresponding  saturation  parameter  fi.  In  figures  5  and  6  the  percent 
difference  and  saturation  parameter  are  plotted  versus  l/\  and  r/ 
respectively.  (In  presenting  the  results,  it  is  more  natural  to  use  the 
normalized  distances,  as  then  the  moments  have  no  explicit  wavelength 
dependence.)  Notice  that  the  percent  difference  goes  to  zero  as  tire 
saturation  parameter  tends  to  one,  as  it  must.  Even  for  low  saturation,  the 
percent  difference  is  still  small.  The  results  are  similar  for  a  vector  von 
Karman  spectrum. 

We  also  find,  numerically,  that  the  estimate  of  is  uncoupled  from  the 
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Normalized  range  iTk  Normalized  range  rtu 


Normalized  range  r//^  Normalized  range  r//. 


Figure  6.  Percent  difference  of  calculated  with  a  nonzero  mean  and  with  a  zero  mean 
versus  r/A  is  in  (a)  and  (b)  for  different  values  of  Z/A  and  Corresponding  saturation 
parameter  is  in  lower  figures  (c)  and  (d).  Calculations  are  in  2D  for  v?  =  0,  d/A  =  0.5, 
n  =  5,  and  a  scalar  von  Kirmdn  spectrum. 


estimates  of  all  the  other  parameters:  r,  I,  and  SNR  (or  NSR).  The 
estimates  of  (p  and  r  will  decouple  provided  that  the  origin  is  taken  at  the 
center  of  the  array  and  that  the  array  is  symmetric;  again  a  linear  array 
with  mirror  symmetry  about  the  origin  is  sufficient.  There  are  no 
conditions  for  the  estimate  of  ip  to  decouple  from  the  estimates  of  the 
remaining  parameters.  If  we  consider  a  simple  two-element  array,  we  can 
readily  see  why  this  occurs. 

Suppose  that  ©a  =  p.  Inspection  of  equation  (D-9)  reveals  that  the 
estimate  of  p  will  decouple  from  the  estimate  of  any  parameter  of  which 
<f>i  is  independent.  Therefore,  as  (l)i  is  independent  of  I,  and  SNR,  the 
estimates  of  these  parameters  decouple  from  the  estimate  of  p.  For  0^  =  r, 
we  have 


_  2M‘^p?k  /  d^2  \ 

a  +  /3  \dp  dp  ) 


o  +  /? 


-  {x[  -P  X2)  sin  p+{y'i-\-  j/g)  cos  p]  . 


(69) 
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Therefore,  the  estimates  of  r  and  ip  will  decouple  provided  x'-^  =  and 
y'l  —  ~2/2- 

6.2  Angle-of-Amval  Estimates  in  Three  Dimensions 

As  we  wish  to  meet  the  conditions  of  the  small-angle  approximation,  we 
take  the  y2;-plane  to  be  the  array  plane  and  the  rc-axis  to  be  near  the 
propagation  axis.  As  noted  earUer,  the  estimates  of  9  and  p  will  be  coupled 
for  this  choice.  As  with  the  two-dimensional  case,  we  find  (numerically) 
that  for  the  estimates  of  9  and  p  to  decouple  from  the  estimates  of  r  and  x, 
we  must  take  the  origin  to  be  at  the  array  center  and  use  a  symmetric 
array  configuration  such  as  a  rectangular  grid  with  mirror  symmetry  in  y 
and  z  about  the  origin.  The  estimates  of  9  and  p  are  found  (numerically)  to 
decouple  from  the  estimates  of  all  the  remaining  parameters.  The  results 
are  independent  of  the  value  of  %. 
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7.  Results 


7.1  Two-Dimensional  Analysis 

It  was  shown  in  section  6.1  that  the  percent  difference  in  the  CRLB  of  ^ 
when  using  a  non-zero  mean  versus  a  zero  mean  is  small.  While  the  CRLB 
of  ip  was  examined  in  detail  in  Wilson  [2],  we  present  the  results  here, 
nonetheless,  as  we  are  using  a  modified  version  of  the  2D  correlation 
functions  fs  and  fy.  The  results  presented  are  for  a  linear  array  with  5 
sensors  spaced  in  equal  intervals  of  d.  The  origin  is  taken  to  be  the  array 
center  so  that  the  estimate  of  ip  decouples  from  the  estimates  of  all  other 
parameters.  As  the  CRLB  for  M  independent  and  identically  distributed 
datasets  is  1/ \/M  times  the  CRLB  for  one  dataset,  all  results  are  presented 
for  M  =  1. 

The  CRLB  of  93  for  y)  =  0  is  plotted  versus  and  I /A  in  figures  7  and  8  for 
SNRs  of  10  and  20  dB,  respectively,  and  for  fixed  values  of  the  other 


lO-**  10^ 


Figure  7.  2D  analysis  of  versus  and  l/X  for  v?  =  0,  r/A  =  500,  d/X  =  0.5,  n  =  5, 
SNR  =  10  dB,  and  a  scalar  von  K4rm^n  spectrum. 
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Figure  8.  2D  analysis  of  versus  and  Z/A  for  =  0,  r/A  =  500,  djX  —  0.5,  n  —  5, 
SNR  =  20  dB,  and  a  scalar  von  Karman  spectrum. 

parameters.  A  scalar  von  Karman  spectrum  is  used.  A  peak  is  evident  in 
at  the  larger  values  of  and  smaller  values  of  //A.  In  this  region  both 
exp  [-OL{pij)r\  and  exp  (-27r)  are  small.  In  fact,  in  the  limit  that 
exp  [-alpij)r]  and  exp  l-2'yr)  simultaneously  approach  zero,  ->  cx)  as 
C  — >  cr^In  and  /i  — >  0.  For  small  values  of  and  oil/X,  both 
exp  \-oi{pij)r\  and  exp  {-2'yr)  are  approaching  the  maximum  value  of  1.0, 
and  hence  is  approaching  the  limit  of  {k'^cP  cos^  ppl/a^)  for  no 
turbulence.  The  value  of  increases  with  decreasing  SNR. 

In  figures  9  and  10,  (Xip  for  (p  =  0  is  plotted  versus  d/X  and  again  for  two 
values  of  the  SNR,  fixed  values  of  all  the  other  parameters,  and  a  scalar 
von  Karmen  spectrum.  In  figure  11,  for  =  0  is  plotted  versus  d/ X  and 
r/A.  For  all  cases,  the  peak  in  occurs  when  botii  exp  [-ot{pij)r\  and 
exp  {—2'yr)  are  small. 

In  figure  12  we  examine  (T(^  at  (/?  =  0  as  a  function  of  the  normalized 
propagation  distance  r /A  for  a  couple  of  values  each  of  the  SNR,  c^,  arid 
I  /A.  For  =  10~®,  the  performance  degrades  with  increasing  r/A, 

especially  for  r/A  >  10^.  However,  for  =  10"^,  the  performance 
degrades  almost  exponentially  for  r/A  >  10^  when  l/X  =  100  and  for 
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Figure  9.  2D  analysis  of  cr^  versus  d/X  and  for  y?  =  0,  r/A  =  500,  //A  =  10,  n  =  5, 
SNR  =  10  dB,  and  a  scalar  von  Kdrmdn  spectrum. 


^/A  ^  10^  when  l/X  —  10.  Intuitively,  we  expect  that  for  regions  of  high 
turbulence  the  best  accuracy  will  be  attained  for  short  propagation 
distances. 

The  CRLB  of  for  other  values  of  has  the  same  behavior  as  for  ip  —  Q, 
with  increasing  slightly.  In  figure  13,  is  plotted  versus  ip  for  two 
values  of  each  r/A,  l/X,  and  SNR.  A  scalar  von  Karman  turbulence 
spectrum  is  used.  We  see  that  for  =  lO"^  depicted  in  the  upper  graph 
(a),  cr^  is  most  sensitive  to  the  SNR;  whereas,  for  10”^  depicted  in  (b), 
is  more  sensitive  to  changes  in  the  normalized  length  scale. 

The  results  for  a  vector  von  Karmdn  spectrum  are  similar  to  those  of  the 
scalar  spectrum:  the  shapes  are  nearly  identical,  but  the  value  of  is  for 
the  most  part  larger.  In  figures  14  and  15  the  ratio  is  plotted  versus 
l/X  and  for  SNRs  of  10  and  20  dB,  respectively.  The  superscripts  v  and  s 
respectively  refer  to  a  vector  and  scalar  spectrum.  Note  that  the  largest 
difference  occurs  when  exp  [—OLy{pij)T\  is  appreciably  smaller  than 
exp  [~as{pij)rY 


(deg) 


Figure  12.  2D  analysis  of  cr^  versus  r/A  for  different  values  of  SNR,  and  l/X.  All  curves 
are  for  at  y?  =  0,  d/A  =  0.5,  n  =  5,  and  a  scalar  von  K^rmdn  spectrum. 


<P  (deg) 


Figure  13.  2D  analysis  of  versus  (p.  A  scalar  von  K^rm^  spectrum  is  used. 


31 


Figure  14.  2D  analysis  of  ratio  crj/crj,  versus  l/X  and  for  v?  =  0,  r/A  =  500,  d/\  =  0.5, 
n  =  5,  and  SNR  =  10  dB. 


Figure  15.  2D  analysis  of  ratio  versus  l/\  and  for  y-  =  0,  r/A  =  500,  d/A  =  0.5, 

n  =  5,  and  SNR  =  20  dB. 
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7.2  Three-Dimensional  Analysis 

We  choose  the  a;-axis  to  be  the  near-propagation  axis  and  take  the  sensor 
array  to  be  in  the  y^-plane.  As  noted  earlier,  this  will  result  in  the  coupling 
of  the  estimates  of  tp  and  6.  The  array  geometry  considered  for  this 
analysis  is  a  4  x  4  square  grid  with  spacing  of  d.  The  origin  is  taken  to  be 
the  array  center.  In  this  way,  the  estimates  of  p  and  6  will  decouple  from 
the  estimates  of  r  and  x,  as  well  as  the  estimates  of  the  other  parameters. 
Again  all  results  presented  are  for  M  =  1. 

The  values  of  a^p,  o-q,  and  are  the  same  at  y;  =  0  and  6  =  7r/2  due  to 
symmetry.  In  figures  16  and  17,  cr^  at  99  =  0  and  0  =  7r/2  is  plotted  versus 
l/X  and  for  SNRs  of  10  and  20  dB,  respectively,  and  fixed  values  of  the 
other  parameters.  A  scalar  von  Karman  turbulence  spectrum  is  used.  We 
see  the  same  behavior  as  for  the  two-dimensional  case.  The  peak  in 
corresponds  to  the  region  where  both  exp  [—a{pij)r\  and  exp  {—2'yr)  are 
small.  The  shape  is  similar  to  that  in  2D,  but  the  overall  values  are  smaller, 
especially  for  large  values  of  and  small  values  of  Z/A.  This  is  expected  as 
a  4  X  4  array  in  3D  should  provide  more  "information"  than  a  5-point 


Figure  16.  3D  analysis  of  cr,^  versus  l/\  and  for  ip  =  Q,  0  =  7r/2,  r/A  =  500,  d/X  =  0.5, 
Uy  =  Uz  =  4,  SNR  =  10  dB,  and  a  scalar  von  K4rm4n  spectrum.  At  =  0  and  0  —  7r/2, 

(Tip  ~  Gq  ~  ~  ^^0* 
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10^ 


Figure  17.  3D  analysis  of  versus  IjX  and  for  —  0,  ^  =  1^12,  r/\  —  500,  d/X  =  0.5, 
riy  —  riz  =  4,  SNR  —  20  dB,  and  a  scalar  von  Kdrm^n  spectrum. 


linear  array  in  2D.  The  behavior  for  other  values  of  and  6  is  similar. 

In  figture  18  we  examine  a^at(p  =  0  cind  9  =  7r/2  as  a  fimction  of  the 
normalized  propagation  distance  r/ A  for  a  couple  of  values  each  of  the 
SNR,  and  l/X.  As  with  the  2D  case,  we  see  that  the  performance 
degrades  most  rapidly  with  increasing  r/A  for  =  10"^;  for  the  3D  case, 
however,  the  near  exponential  rise  in  now  occurs  at  slightly  larger 
values  of  r/A. 

As  in  the  2D  case,  the  behavior  of  cr,^  and  ae  for  a  vector  von  Karm^ 
spectrum  is  similar  to  that  of  a  scalar  von  Karman  spectrum:  the  CRLBs  for 
the  vector  spectrum  increase  most  significantly  over  those  for  the  scalar 
spectrum  in  regions  where  exp  [-ay{pij)r]  is  appreciably  smaller  than 
exp  [— as(pij)r].  The  ratio  at  y;  =  0  and  6  —  7r/2  is  plotted  versus  l/X 

and  for  SNRs  of  10  and  20  dB  in  figures  19  and  20,  respectively. 

Figure  21  examines  a^p  and  ae>  as  fimctions  of  and  O',  as  well  as  the 
coupling  between  the  estimates  of  ip  and  9',  for  =  10~®  and  l/X  =  10.  A 
scalar  von  K^rmdn  spectrum  is  used.  From  figure  21(a),  a  contour  graph  of 
(T^  for  9'  versus  ip,  we  see  that  for  small  angles  9')  =  Op{9',  p).  This 
symmetry  was  also  seen  for  small  angles  in  the  no-turbulence  case. 
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Figure  18.  3D  analysis  of  versus  r/A  for  =  0,  0  —  7r/2,  djX  —  0.5,  riy  =  riz  =  4,  and  a 
scalar  von  Kerman  spectrum. 


Figure  21(b)  gives  cr^p  as  a  function  of  (p  for  fixed  values  of  9'.  Now,  ue',  on 
the  other  hand,  is  dependent  upon  6'  but  is  independent  of  p,  as  indicated 
in  figure  21(c).  Again  this  is  the  same  behavior  as  in  the  no-turbulence 
case.  As  the  estimate  of  neither  p  nor  9'  is  coupled  to  the  estimate  of  any 
other  parameter,  we  may  use  equation  (8)  to  gauge  how  strong  the 
coupling  between  the  estimates  of  p  and  9'  is.  We  first  note  that  and 
are  both  dependent  upon  p  and  9'.  In  figure  21(d),  C  is  plotted  versus 
p  for  fixed  values  of  9'.  We  see  that  the  coupling  is  in  fact  weak.  In 
figure  22(a-d)  the  same  is  presented  but  for  —  10“^  and  //A  =  10. 
Figures  21(d)  and  22(d)  appear  identical.  In  fact  we  find,  to  within  the 
numerical  accuracy  of  this  calculation,  that  is  the  same  as  for  the  no 
turbulence  case  (see  equations  (57)  and  (60)),  regardless  of  the  the  value  of 
r/A,  d/X,  SNR,  I  /A,  or  and  regardless  of  the  von  Karman  spectrum 
used. 
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10"^  10*'  //A 


Figure  19.  3D  analysis  of  ratio  versus  //A  and  for  v?  =  0,  0  =  7r/2,  r/A  =  500, 

d/X  =  0.5,  Uy  =  riz  =  4,  and  SNR  —  10  dB. 


Figure  20.  3D  analysis  of  ratio  (tJ/ct*  versus  //A  and  for  v?  =  0,  0  =  7r/2,  r/A  =  500, 
d/X  =  0.5,  Uy  —  riz  ~  4,  and  SNR  =  20  dB. 


SQRTCRLB  of  elevation  (deg)  Elevation  (deg) 


(a)  SQRTCRLB  of  azimuth  (deg)  (b)  SQRTCRLB  of  azimuth 


(c)  SQRTCRLB  of  elevation 


X  10'^  Coupling  coefficient 


Figure  21.  (a)  Contour  of  for  6'  versus  (p.  (b)  Plot  of  cr^  versus  p  for  fixed  6'.  (c)  Plot  of 
versus  O',  (d)  Plot  of  C  versus  p  for  fixed  O'.  All  curves  are  for  r/X  ~  500,  l/X  =  10, 
=  10"^,  d/X  =  0.5,  ny  —  riz—  4,  SNR  =  10  dB,  and  a  scalar  von  Karman  spectrum. 


(a)  SQRTCRLB  of  azimuth  (deg) 


(b)  SORTCRLB  of  azimuth 


Azimuth  (deg) 


Azimuth  (de^) 


(c)  SQRTCRLB  of  elevation 


(d)  Couplins  coefficient 


Independent  of  azimuth 


^  4.95 


5  10 

Elevation  (deg) 


5  10 

Azimuth  (deg) 


Figure  22.  (a)  Contour  of  for  &  versus  v?.  (b)  Plot  of  versus  ^  for  fixed  (c)  Plot  of 
oq^  versus  0'.  (d)  Plot  of  C  versus  v?  for  fixed  All  curves  are  for  r/A  =  500,  //A  =  10, 
^2  =  10  ^  ,dl\  —  0,5,  riy  —  Uz  —  4,  SNR  =10  dB,  and  a  scalar  von  K^rm^n  spectrum. 


7.3  Cramer-Rao  Lower  Bound  of  Other  Parameters 

Though  we  have  not  examined  the  CRLB  of  the  other  parameters  in  this 
report,  a  brief  comment  is  in  order.  For  large  values  of  the 
index-of-refraction  variance  and  small  values  of  the  normalized  length 
scale,  we  find  that  for  0  G  [Z/A,  r/A] ,  <Je/0  may  be  less  than  one  for  a 

sufficient  number  M  of  independent  and  identically  distributed  datasets. 
However,  for  other  values  of  and  Z/A,  uq/Q  >  1.  Therefore,  the  regions 
in  which  cr^  is  strongly  dependent  upon  0  correspond  to  the  regions  in 
which  a@ /0  is  smallest;  similarly  the  regions  in  which  is  nearly 

independent  of  0,  correspond  to  the  regions  in  which  0  cannot  be 
estimated  from  the  model  presented  here.  Due  to  the  complexity  of  the 
analysis  of  ere,  it  is  left  for  another  date. 
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8.  Summary  and  Conclusions 


We  have  examined  the  performance  boimds  of  acoustic  arrays  operating 
in  atmospheric  turbulence  with  fluctuations  described  by  a  von  K^rm^n 
spectrum.  The  analysis  featured  three  main  improvements  upon  earlier 
work: 

•  The  performance  bounds  were  generalized  to  weak  as  well  as  strong 
scattering. 

•  Multiple  unknowns,  such  as  the  propagation  distance  of  the 
wavefront  and  the  turbulence  parameters,  were  incorporated. 

•  AOA  estimates  for  three-dimensional  problems  (i.e.,  two  bearing 
angles)  were  considered. 

Our  primary  interest  in  this  report  was  to  analyze  the  CRLBs  of  the  AO  As. 
For  an  incident  plane  wave,  we  have  found  that  an  appropriate  choice  of 
coordinate  system  and  array  geometry  leads  to  the  decoupling  of  the 
estimates  of  the  AO  As  from  the  estimates  of  the  other  parameters:  the 
normalized  propagation  distance  (r/A),  SNR,  turbulence  parameters  (Z/A 
and  c^),  and  phase  angle  (x).  In  order  to  remain  consistent  with  the 
small-angle  approximation,  we  had  to  choose  a  coordinate  system  that 
resulted  in  the  coupling  of  the  estimates  of  the  azimuth  and  zenith;  the 
coupling,  however,  was  small.  We  have  found  that  the  CRLBs  of  the  AOAs 
increase  significantly  for  large  values  of  the  normalized  propagation 
distance  and  index-of-refraction  variance.  However,  for  small  values  of 
the  index-of-refraction  variance  and  normalized  propagation  distance,  the 
SNR  is  the  limiting  factor. 

The  results  in  this  report  demonstrate  that  scattering  by  atmospheric 
turbulence  sigmficantly  affects  the  performance  of  acoustic  sensor  arrays. 
Future  efforts  should  attempt  to  incorporate  the  additional  phenomena  of 
ground  reflections  and  refraction  by  atmospheric  wind  and  temperature 
gradients.  These  phenomena  will  likely  have  considerable  impact  on  the 
ability  to  estimate  the  elevation  angle.  However,  numerical  techniques 
will  probably  be  required  to  model  these  complications. 
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Appendix  A.  Nomenclature 


A.l  Symbols 

£  Complex 

=1=  Complex  conjugate 

=  Defined  as 

e  Element  of 

(•)  Ensemble  average  or  expectation  value 

t  Hermitian  adjoint  (complex  conjugate  transpose) 

In  Identity  matrix,  n  x  n 

9  Imaginary 

~  On  the  order  of 

3?  Real 

T  Transpose 


A.2  Abbreviations 

AOA 

AOB 

CRLB 

FI 

LOS 

MLE 

MSE 

PDF 

SQRTCRLB 

2D 

3D 


Angle  of  arrival 
Angle  of  bearing 
Cramer-Rao  lower  boimd 
Fisher  information 
Line  of  sight 

Maximum  likelihood  estimator 

Mean-square  error 

Probability  density  function 

Square  root  of  the  Cramer-Rao  lower  boimd 

Two  dimensions  or  two-dimensional 

Three  dimensions  or  three-dimensional 


Appendix  B.  Derivatives  of  Covariance  Matrix  and  Mean 


In  order  to  calculate  the  FI  matrix  (equation  (5)),  we  need  the  derivatives 
of  the  covariance  matrix  C  and  the  mean  vector  with  respect  to  the 
unknown  parameters  ©  =  [<p  0  r  Z  x]  The  dependence  of  C  and  fj, 
on  these  parameters  is  given  in  section  4.1.  The  derivatives  of  the 
components  of  the  covariance  matrix  (Cu  and  Cij,  i  /  j)  and  mean  vector 
(/ij)  with  respect  to  these  parameters  are  given  in  the  following 
subsections. 


B.l  Unknovm  Parameter  ip 

The  derivatives  with  respect  to  the  azimuth  p  are 

dCii 
dtp 
dCjj 
dp 
dpi 
dp 

where 

o  / 

-^-1  =  k  sin  p  sin  0  +  yl  cos  p  sin  0) 

dcpij  _  d^  _  d(t)j 
dp  dp  dp 


=  0 


-  i  c-  i 

.  d(t)i 


B.2  Unknown  Parameter  6 

The  derivatives  with  respect  to  the  zenith  (decimation)  0  are 


dCii 

d0 

dCi 


=  0 


d0 


=  1 


dpi 

d0 


.  d(j)i 


3 


(B-1) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 


(B-6) 

(B-7) 

(B-8) 
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where 


=  -k{x[  cos  ip  cos  6  +  y\  sin  cos  0  -  z\  sin  9) 
d(l)ij  _  d4>i  d(})j 

'W 


(B-10) 


B.3  Unknown  Parameter ; 


The  derivatives  with  respect  to  the  propagation  distance  r  are 


J^ii  _L  ; 


+  i  Cij  ,  i^j 


i*,  . 


dr  dr 


+1^  Hi, 


(B-11) 


(B-12) 


(B-13) 


where 


-a(oo) 


(B-14) 


=  -7M 


(B-15) 


B.4  Unknown  Parameter  for  Noise 


(B-16) 

(B-17) 


We  can  treat  either  the  noise  variance  (ct^  )  or  the  signal-to-noise  ratio 
(SNR)  as  the  unknown.  First  let  us  suppose  that  is  the  unknown.  Then 


(B-18) 


{B-19) 


(B-20) 
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Alternatively,  we  can  take  the  SMR  as  the  unknown.  Recall  that 
SNR  =  -  Let  0  be  the  parameter  symbol  for  the  SNR.  Then 


dCii  dal 


de 


dCi, 

de 

dp 


de 

=  0,  i^j 


^=0, 


where 


_  PO  _ 

de  02  0  ■ 


(B-21) 

(B-22) 

(B-23) 


(B-24) 


B.5  Unknown  Turbulence  Parameters 

Consider  the  derivatives  with  respect  to  the  turbulence  parameters  I  (the 
length  scale)  and  (the  index-of-re£raction  variance).  Let  0  G  [/ ,  be 
the  unknown  parameter.  Then 


dCii  dp? 


de 

dCij 

de 


de 

2  dajpij)  a(pij)  r  _  ^ 

ae  de 


^  ^  J<t>i 

ae  ae 


(B-25) 

(B-26) 

(B-27) 


where 


ap^ 

ae 

dp 


„  da{oo)  .2 
“5 - —  P 


ae 

r  aa(oo) 

2  ae 


p. 


(B-28) 

(B-29) 


Recall  that  a{p)  =  2'Kk‘^  [/(O)  -  /(p)],  where  /  is  the  2D  correlation 
function  for  either  a  scalar  or  vector  von  Karman  spectrum  as  given  in 
equations  (18)  and  (19),  respectively.  The  derivatives  of  /«  and  with 
respect  to  I  and  are  given  in  the  following  appendix. 
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Appendix  C.  Derivatives  of  /(p)  with  respect  to  and  I 


The  scalar  arid  vector  2D  correlation  fimctions,  /«  and  /y,  are  given, 
respectively,  in  equations  (18)  and  (19).  In  this  section  we  calculate  the 
derivatives  of  fs  and  /„  with  respect  to  the  index-of-refraction  variance 
and  with  respect  to  the  turbulence  length  scale  1.  The  limits  of  /  (p),  for 
p  ^  0  and  p—^oo,  are  also  calculated  as  well  as  their  derivatives. 


C.1  Derivative  of  f{p)  with  respect  to 

The  derivative  of  /,  where  /  is  either  /«  or  ,  with  respect  to  the 
index-of-refraction  variance  is 


Q(-2  ^2 


(C-30) 


C.2  Derivative  of  /(p)  with  respect  to  p  7^  0 

The  derivative  of  the  scalar  2D  correlation  function  /s(p),  for  p  7^  0,  with 
respect  to  the  turbulence  length  scale  I  is  given  by 


dfs 

dl 


dl 


\/^r(i/3) 


v^r(i/3) 

2?2 

v^r(i/3) 

2c2 

V^r(i/3) 


(t)  +  21  (7)  ( 


2^2  /p\5/6 

v/^r(i/3)  \2i) 


[Kn/e  (y)  +  y^i/e  (y)_  • 


(C-31) 


In  going  from  the  second  line  to  the  third  line  and  from  the  forth  line  to 
the  fifth,  the  recxursion  relationships 

K,-i{x)  +  K^+i{x)  =  -2KUx)  (C-32) 

Ot/ 

K,-i{x)  -  (a:)  =  -—KA^)  (C-33) 
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and  the  relationship 

K^{x)  =  K^uix)  (C-34) 

have  been  used,  (see,  for  example,  Arfken  [13]  or  Gradshteyn  and  Ryzhik 
[14]).  The  intermediate  algebraic  steps  are  not  shown. 

The  derivative  of  the  vector  2D  correlation  function  fy{p),  for  p  ^  0,  with 
respect  to  the  turbulence  length  scale  I  may  be  written  in  the  form 


^  ^  I  [,-5/6 fP\] 

di  di  0?r(i/3)V2y  dii  '/nJj 

V^r(i/3)  (2)  “6 (?)  “ (?)] 

v/)Fr(i/3)  (S)  {~6  (?)  ^  21  (?)  +  (? 


dl  v/^r(l/3)  V2//  \  6  21 

Ar(i/3)  (i)  {6^"/®(?)'^l  .f^'/®(?)+^iV6(?) 


vAFr(l/3) 


~(l)  [^^/®(?)+^^/e(?)]} 


(C-35) 


Again  the  recursion  relationships  have  been  used  in  going  from  the 
second  line  to  the  third  and  from  the  forth  line  to  the  fifth.  The 
intermediate  algebraic  steps  are  not  shown. 

C.3  Evaluation  of  /(O)  and  Derivative  of  /(O)  with  respect  to  I 

We  wish  to  evaluate  fs{p)  and  fy{p)  in  the  limit  that  p  ^  0.  Consider 


/„(0)  =  lira  C  ([)  '  (^)  -  lK,/e 


(C-36) 

(C-37) 


where 


2V60i:r(l/3)  ■ 

For  X  E  5R  and  v  >  0,  (see,  e.g.,  Arfken  [13]), 

,.  ^  ,  r(z/)2"-i 
urn  A„(x)  =  — - 

i_o  '  x‘' 


(C-38) 


(C-39) 
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Therefore,  for  a;  €  3?  and  u  >  0, 


oo  if  A  <  1/ , 

lim  =  <  2‘'-'^T{u)  if  X  =  u,  (C-40) 

,0  if  A  >  1/ . 


It  follows  that 


/,(0)  =  C2-v«r(5^  = 


%AFr(i)  ;3(1,1) 


fv{0)  ) 


where  (3  is  the  beta  function,  I3{x,  y)  =  r(x)r(i/)/r(a;  +  y).  Thus 


dm  ^  _  /(O) 

dl  i 


(C-41) 


(C-42) 


C.4  Evaluation  of  /(oo) 

We  wish  to  evaluate  fs(p)  and  fv{p)  in  the  limit  p  oo.  We  note  that  in  the 
asymptotic  expansion  of  Kt,{x),  x  G  SR,  the  leading  term  is  proportional  to 
e“® /v^.  As  such  the  exponential  term  will  dominate  in  fs  and  /„.  Hence 

fs{oo)  =  fv{oo)  =  0 .  (C-43) 

We  refer  the  reader  to  [13]  and  [14]  for  comments  regarding  the 
asymptotic  expansion  of  Ku{x). 


Appendix  D.  Simple  Two-Element  Array  Model 


Consider  the  simple  case  of  a  two-element  array,  and  let  ©>  and  0j,  be 
unknown  parameters.  Let  us  write  the  covariance  matrix  as 


C  = 


a  b 
b*  a 


(D-1) 


where  b  =  and  a,  /?  G  Ji.  Let  us  write  the  mean  vector  as 


Ml 

LM2j 


=  M 


e 


(D-2) 


tr  C 


It  follows  that 

-1  dc  ac  \  2 


dOx  dej  (detC) 


3? 


Y  da  i,*  db  \  f  da  db  \ 

^  dQx)  V  50^  ^  dQu) 

da 


db 

•h  1  a-;:^^ - b 


db*  _  da 
a^r^ - b 


Y0a  ''dexJVde^  “50^ 

/  2  /i2\  / 


(D-3) 


'  ^\a0Aa0^  a0Aa0„^ 

^  da  dd  da  dp  \  2  /  2  /o2'|  d<t>\2  d(t>i2 

P-4) 


and 


3fi 


d/x  dfi 


dQx 


d9^ 


=  sf?  /  ^  _ 

\detC  [a0A  V  90^ 


.  ^M2\  ^ 

dej  dGx  V  90: 


9fl2  ,, 
- 0 


9fii 


a?  —  Y 


a0A  a0j. 


50^ 

9<f>2  9(j)2 


(D-5) 


/  d(pi  9(j)i 

^  V^0A  9Qu  9Gx  90 


501  502  501  502 

50a  50^  50^  50a 


I)]}- 


(D-6) 


Recall  for  a  Gaussian  process,  the  FI  is  given  by  equation  (5).  Thus,  for 
^  =  A  or  I/, 
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For  the  special  case  when  da/dQ\  =  dP/dQx  =  dfi/dQx  =  0,  but  there  are 
no  restrictions  on  the  derivatives  with  respect  to  0^, 


_  2M/32  /dduV  2Mij?  [ 
a?~l3^  \dex)  ^  a2-/?2 

^  2Mp^  9012  9012  2Mfi? 

~  a?  -  /32  dQx  dQ^  ^  a2  - /52  [“ 


-2P 


d(j)i  dd)2  1 

5eA^J 


/  d(t>i  d(j)-i  d(l>2  d<j>2  \  /  501  502 

\dex  50,  50a  50  J  “  ^  [dQx  50, 


(D-8) 

501  502  y 

50,  ’ 

(D-9) 


and  J„  is  given  by  equation  (D-7)  with  ^ 
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